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The single electron Green's function of the one-dimensional Tomonaga-Luttinger model in the pres- 
ence of open boundaries is calculated with bosonization methods. We show that the critical ex- 
ponents of the local spectral density and of the momentum distribution change in the presence 
of a boundary. The well understood universal bulk behavior always crosses over to a boundary 
dominated regime for small energies or small momenta. We show this crossover explicitly for the 
large-U Hubbard model in the low-temperature limit. Consequences for photoemission experiments 
are discussed. 
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There is currently great interest in "Luttinger liquid 
physics" sparked by a new generation of exper- 

iments on low-dimensional electron structures. Exam- 
ples include measurements of the point contact tunneling 
conductance between two fractional quantum Hall edge 
channels |3j and high-resolution photoemission studies of 
quasi one-dimensional metals |^,|^. Additional interest 
stems from the fact that the Luttinger liquid - i.e. the 
low-energy, long-wavelength physics of interacting elec- 
trons in one dimension (ID) - provides us with the only 
known non-Fermi liquid phase with unbroken symmetry. 
For this reason the notion of a "Luttinger liquid" (LL) 
has played a prominent role in studies of generic features 
of correlated electron systems, and it has been suggested 
that some of its properties (e.g. anomalous propagators, 
spin-charge separation) carry over to higher dimensions 

So far little attention has been paid to the effect of 
boundaries on the spectral properties of Luttinger liq- 
uids (except in the context of spin-chains |Q). This is 
surprising since boundary effects are bound to be impor- 
tant in several of the proposed laboratory realizations. 
For example, in a recent series of experiments on the 
Bechgaard salts (TMTSF)2X (X is a counter-ion), the 
Hel and II photoemission spectra were measured in the 
metallic phase Q . This class of materials are known to 
exhibit strong electron correlations and are prime candi- 
dates for LL-behavior. The experimental results suggest 
an anomalous suppression of the spectral weight close to 
the Fermi level. It was argued in ||] that consistency with 
the LL-scenario requires long-range electron-electron in- 
teractions, but as we will show here, boundary effects 
also deplete the spectral weight and may influence the 
observed data. 

Specifically, we have examined the effect of an open 
boundary on the local single-particle spectral density of 
a spinful LL. We have also studied the momentum dis- 
tribution in a finite system with open boundary condi- 



tions. In both cases we observe a strong influence from 
the boundaries, which cause novel scaling behavior with 
energies and momenta close to the Fermi level. 

As our model we take an extended version of the 
Tomonaga-Luttinger (TL) Hamiltonian, describing the 
low-energy limit of locally interacting ID electrons 
defined by the Hamiltonian density: 
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are the chiral Fermion cur- 
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rents of the left and right moving components of the elec- 
tron field ^'ct(x), expanded about the Fermi points ±fci?: 
^^{x) = e-*'=^^'(/'L,rT(2:)-he^'=^^Vii,<T(a:). The Hamihonian 
(|l|) describes left- and right moving relativistic Fermions 
in (l-l-l)-dimensions which interact via forward scatter- 
ing without spin flip (51,52, and 53) or with spin fiip 
(54). "Umklapp" processes are suppressed away from 
half-band-filling, so this model provides a complete pic- 
ture of possible local interactions in this case. (Trivial 
forward scattering terms which can be absorbed by re- 
defining the Fermi velocity have not been explicitly writ- 
ten out.) 

The TL Hamiltonian is conveniently bosonized |^ by 
introducing charge and spin currents and the correspond- 
ing bosons 4>c and (f)s with conjugate momenta H^ and H^, 
respectively: 
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and accordingly for right-movers. The resulting theory 
describes separate spin and charge excitations moving 
with velocities Vc = ^ + ^ and Vs = ^ — respec- 
tively: 
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+ 54 const, cos ^/^TTips, (3) 



where 5c = 5i + 92 and gs —92 — 91- The charge inter- 
action 5c can be absorbed into the free Hamiltonian by a 
simple rescahng of the charge boson, but the spin inter- 
actions 5s and 54 obey KosterUtz-Thouless renormahza- 
tion group equations |l0[ with flow hues along hyperbolas 
5s —54 = const, (to lowest order). For gs > — I54I the spin 
sector develops a gap in the low energy, long wave-length 
limit, but for gs < — I54I the system flows to a stable fixed 
point g* = -^^g'^s - 9h 9% = 0. For 5^ = -I54I the in- 
teraction corresponds to one single marginally irrelevant 
operator, so that gt — 9t — 0. If the flow to a stable 
fixed point occurs, we can rescale the bosons by a canon- 
ical transformation to obtain a free theory [v = s, c) 

(P.^K.cj,,, ^n,/x„ (4) 

where to first order in the coupling constants 

= l~g*jATrVs, Kl^l-gjAiiv,. (5) 

We now consider a semi-infinite system with an open 
boundary condition at the origin and thus require the 
electron field '^a{x) to vanish at x — 0. This im- 
plies tl>L,a{^) = — ^_R,cr(0), or in terms of the bosons 
H..M = -(t>R,M + W/2i^c and <^L,s(0) = -0fl,s(O), 
which allows an analytic continuation of the left-movers 
onto the negative half axis in terms of right-movers 



<l>L,u{x, t) = -~4>R.u{—X: t) + const., X < Q 



(6) 



(i^ = s,c, const. = 0, -y7r/2_ft'c)- We can therefore 
describe the theory in terms of left-movers only which 
live on the full complex plane without an explicit bound- 
ary condition. Using this formalism, the single electron 
Green's function can be calculated in a straightforward 
way: 

<^Ux,t)^^iy,0)> 

- e''"'^'=+y'>G{x,~y,t) - e~'''^^''+y'>G{-x,y,t), (7) 

where the chiral Green's function G{x,y,t) is a product 
of spin and charge contributions 
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Here, x and y denote the distance from the boundary 
(a; = 0) and the time carries an implicit ultraviolet cut- 
ofl[ t — ie. We can see that in the limit xy \{x — v)^ — 



v^t I the last factor in equation (Q) goes to unity and we 
recover the known bulk correlation function ||^ (as we do 
in the non-interacting case Kc — = 1). 

To understand the physical implications of the bound- 
ary correlation function we study the local spectral den- 
sity N{u), r), which is given in terms of the Green's func- 
tion as 
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N{cj,r)^- e^-\{^Ur,0),^Ar,t)})dt, (9) 
^~ J-00 

where lu is measured relative to the Fermi energy and r is 
the distance from the boundary. Without the boundary, 
this integral can be done exactly with the result that 
the spectral density scales at the Fermi surface as N{u!) ex. 
Lj°' , where the exponent in the bulk is given by 



a = {Kt 



kI + k;^)/a-i. 
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However, the boundary clearly has an effect on this ex- 
ponent, and simple power counting shows that we expect 
a crossover to a boundary dominated regime for uj < v^/r 
with a novel exponent as = {K^'^+K^'^)/2—l. Interest- 
ingly, the boundary exponent as therefore always dom- 
inates for sufficiently small uj. Moreover, we notice that 
the last two terms in equation (^ make a contribution 
which oscillates at twice the Fermi wave- vector and drops 
off with the distance from the boundary proportional to 
gi2fcj^r^-(_ff^+_fcj/2 rpj^-g contribution is reminiscent of a 
Friedel oscillation, although it can probably not be ob- 
served directly, since experimental measurements of the 
density of states (e.g. photoemission) will average over 
several lattice sites. We therefore ignore those "Friedel" 
terms in the following calculations. 

As an example, we consider the low temperature Hub- 
bard model away from half-filling which is well under- 
stood in terms of the TL model In this case, the 
SU(2)-invariance forces Ks — 1, and it is known from 
Bcthc ansatz calculations that ^ 1/2 as C/ ^ 00 
]l2| . From these numbers, the well known result N(lu) cx 
|[j|i/8 for small uj follows immediately for the bulk regime 
LiJ 3> Vc/r. In the presence of the boundary, however, 
we crossover to the boundary exponent as — 1/2 for 
UJ < Vc/r. After rescaling the variable of integration 
in equation (^, we see that the spectral density is a 
function of the scaling variable ruj only (up to an over- 
all constant): N{r,Lj) = r~^/* fir^)- After a deforma- 
tion of the integration contour, a numerical integration 
of equation (||) is straightforward (we subtract the di- 
vergent part). The results in figure (|l]) clearly show the 
crossover from boundary behavior for ruj/vc < 1 with 
exponent as = 1/2 to bulk behavior for ruj/vc > 1 with 
exponent a — 1/8 (the corresponding power-laws are su- 
perimposed). The observed oscillations in figure (|^) are 
an interesting secondary effect, which vanish asymptot- 
ically as sin(2a;r/uc)(w?')~^^/^^, (but those are not due 
to the "Friedel" terms in equation (^ which have been 
neglected) . 
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FIG. 1. The spectral density as a function of rcu in arbi- 
trary units. The corresponding power-laws for a — 1/8 and 
Ob ~ 1/2 are also shown. 

Our results have direct relevance for photoemission ex- 
periments on quasi-lD metals. At low temperatures the 
photoemission intensity I{uj) is proportional to the lo- 
cal spectral density N{uj,r), integrated over the region 
of escaping electrons and weighted by the Fermi-Dirac 
distribution 

lico) ^ J dr fFDiu^)N{co,r) . (11) 

(We neglect the small thermal shift in N{lu, r) at low 
temperatures.) 

In a boundary dominated region, I{uj) is seen to be 
dramatically reduced at the Fermi level compared to a 
"bulk region" . Moreover, the finite energy resolution A 
of the photon lines effectively introduces an averaging 
over the "true" spectral density 

J(w)obs = / e-(--")'/2^'/(x) dx. (12) 

This averaging completely wipes out the power-law sin- 
gularities in either the bulk or the boundary case as 
shown in figure (||), where we plotted I{uj)ohs in arbi- 
trary units at T = 50K for boundary {as ~ 1/2) and 
bulk [a — 1/8) regimes, respectively, assuming an exper- 
imental resolution of A = 20meV (experimental values 
according to ||^). The corresponding three-dimensional 
case {a = 0) is also shown for comparison. In the neigh- 
borhood of the Fermi level the observed boundary dom- 
inated spectral density appears to be depleted with an 
exponent of one or larger (compared to the exponent 
a = 1/8 of the bulk spectral function without tempera- 
ture or averaging effects). 
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FIG. 2. The predicted intensity I^hs in arbitrary units as a 
function of uj for boundary and bulk cases (i.e. for power-laws 
with Qs = 1/2 and a = 1/8, respectively). The correspond- 
ing three-dimensional case (a = 0) is also shown. Temper- 
ature (T — 50K) and finite resolution (A — 20meV) effects 
have been taken into account according to equations (O) and 
®- 

In experiments the condition for boundary behavior 
UJ < Vc/r will be fulfilled over an energy range io ~ Ep/£, 
where £ is the distance from the boundary in units of the 
lattice spacing. This means that if the effective length 
i of the ID chains is about one-hundred lattice spac- 
ings, the boundary effects are observed over a region of 

lOOmeV around the Fermi energy (e.g. corresponding 
to an "impurity" density of 1% close to the cleaved sur- 
face) . Turning to the photoemission data reported in [|j , 
a combination of boundary effects and a finite experimen- 
tal resolution go a long way to account for the observed 
suppression of spectral weight at the Fermi level. The 
experiments indeed suggest a scaling /(a;)obs esc uj°'°^'', 
with aobs > 1 extending, however, over a larger energy 
interval, and some additional mechanism (long-range in- 
teractions [^,0 or electron-phonon coupling ||l3|) may 
have to be invoked to fully explain the data. Another 
way of directly examining the boundary region would be 
to probe the single crystal face where the open ends of 
the chains are. Since the electrons have an escape depth 
of only a few lattice spacings, the boundary exponent 
should be observable over a much larger energy range. 

As a second example of the effect of open boundaries, 
we consider the momentum distribution of a finite system 
with length L and open boundary conditions, which can 
be expressed in terms of the chiral Green's function 

1 

n{k,L)=— / dx dy cosk{x — y) G{x,y,0), (13) 
L Jo 

where k is measured relative to the Fermi wave-vector 
hp (we also implicitly subtract n(0,L) to remove any 
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divergences). The bulk behavior can be determined by 
taking L ^ oo and by simple power-counting we see that 
n(fc,oo) oc |fc|", where a is again given in equation (10). 
However, boundary effects will be present, and moreover 
we have to consider that we are now dealing with a finite 
system (i.e. open boundary conditions both at a; = and 
X = L). One way of determining the correct correlation 
functions in this case is to assume conformal invariance 
in the complex plane z' — x + iv^r (which is justified 
in the low energy, long-wavelength limit and with decou- 
pled spin and charge sectors). We then simply perform 
a conformal transformation onto a cylinder with circum- 
ference 2L: z' ^t^rz/L jg^f^gj. ^jjg analytic continuation 
in equation (|^)]. Since we know the chiral Green's func- 
tion in the plane (z^, we immediately obtain the result 
for the finite case |14[| (using the transformation rules for 
primary fields) 



G{x,y,t) cx 



nf 2L . TT{Vi,t 



y) 



' 2L T^iv^t 

— sm 
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X - 
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(14) 



sin ^ sin ^ 



■K{v,^t-\-x-\-y') . ■K{v,,t — x — y) 

. 2L sm ^ 



We expect the critical behavior to be changed dramat- 
ically by this transformation, but now the exponent can- 
not be determined by simply counting the powers in the 
finite Fourier transform. Again, the momentum distribu- 
tion is a function of a scaling variable kL only (up to an 
overall constant): n{k,L) — L^" f{kL). We are mostly 
interested in the apparent critical exponent ctapp near the 
Fermi wave- vector n{k,L) cx |fc|"''pp. This apparent ex- 
ponent will change slowly depending on the scale kL at 
which we probe the system, so it is useful to define a scale 
dependent exponent in terms of the logarithmic deriva- 
tive 

aapp(fci)^^|(fc,L). (15) 

As our example we consider the case of the large- C/ 
Hubbard model again. After calculating ctapp numeri- 
cally as shown in figure (|^), we find that the crossover 
is extremely slow. Even for huge values of kL ~ 10^ we 
are still considerably far away from the accepted bulk 
exponent a — 0.125. Since k is assumed to be small 
compared to kp, this means that even macroscopic sam- 
ples of several centimeters will have a finite-size domi- 
nated momentum distribution near kp- A slow crossover 
is also observed for a periodic (but finite) system, which 
is shown for comparison in figure (^. The observed be- 
havior is therefore mostly a finite size effect, but systems 
with an open boundary condition show an even slower 
crossover. 
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FIG. 3. The apparent exponent tJapp of the momentum 
distribution for a finite system as a function of kL according 
to equation (p^). As kL — » oo, the bulk exponent a = 0.125 
is slowly approached. 

In conclusion we have shown that an open boundary 
has a pronounced influence on the observed critical ex- 
ponents in quasi one- dimensional metals. The spectral 
density was shown to be always dominated by the bound- 
ary exponent for frequencies close to the Fermi energy, 
which has direct consequences for the interpretation of 
photoemission experiments. The momentum distribution 
exhibits strong flnite size and also boundary effects close 
to the Fermi wave- vector and the crossover to "bulk" be- 
havior is extremely slow. 
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